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Abstract 

Using a large sample of D + — > K~K + ir + decays collected by the FOCUS photo- 
production experiment at Fermilab, we present the first non-parametric analysis 
of the K~ir + amplitudes in D + — ► K~K + ir + decay. The technique is similar to 
the technique used for our non-parametric measurements of the D + — > K e + v 
form factors. Although these results are in rough agreement with those of E687, 
we observe a wider S-wave contribution for the K (1430) contribution than the 
standard, PDG [1] Breit-Wigner parameterization. We have some weaker evidence 
for the existence of a new, D-wave component at low values of the K~ir + mass. 



1 Introduction 

This paper describes a non-parametric measurement of the K~n + amplitudes 
present in the decay D + — > K~K + ir + . Charm decay Dalitz plots are tra- 
ditionally fit using the isobar model where the amplitude is represented by 
a sum of known Breit-Wigner resonances multiplied by complex amplitudes 
along with a possible non- resonant term [2] . Often this approach gives a good 
qualitative representation of the observed populations in the Dalitz plot. The 
isobar approach, however, does not automatically incorporate some important 

1 See http://www-focus.fnal.gov/authors.html for additional author information. 
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theoretical constraints. If all final state interactions are dominated by the two- 
body resonant system, with negligible contribution from the third body, then 
two-body unitarity is violated by the isobar model. These potential unitarity 
violations are particularly severe for the case of broad, overlapping resonances. 

An alternative, K-matrix formalism spearheaded by the FOCUS collaboration 
[6] for charm meson decays, is designed to automatically satisfy unitarity but 
it is a difficult analysis that requires considerable input from low energy scat- 
tering experiments on the form and location of K-matrix poles. In particular 
the FOCUS analysis of the D + , — ► 7r + 7r _ 7r + final state used the K-matrix 
description for just the broad dipion states, but the narrow P-wave states were 
incorporated as isobar contributions. 

This work takes a considerably different approach. We directly measure the 
K~7i + spin amplitudes as a function of m^-^ mass by "projecting" them 
based on the decay angular distribution. Our projective weighting technique, 
described in Section 3, is very similar to that used to make non-parametric 
measurements [5] of the q 2 dependence of the helicity amplitudes in D + — > 
K~n + £ + v. In our earlier work, each helicity amplitude is projected by making 
a weighted histogram of q 2 using special weights designed to block all other 
amplitude contributions. The projective weighting technique is therefore an 
intrinsically one dimensional analysis. The D + — > K~K + tt + final state - in 
principle - is influenced by three amplitudes: K'tt + , K~K + , and even K + n + . 
We have chosen the D + — > K~K + tt + final state as a first test case since the 
E687 isobar analysis [3] concluded that the observed D + — > K~ K + n + Dalitz 
plot could be adequately described by just three resonant contributions: 4>tt + , 
K+K , and K+K (1430). Although </>7r + is an important contribution, the 
is a very narrow resonance that can be substantially removed by placing a 
lower cut on the tuk+k- mass (such as rrix+K- > 1050 MeV/c 2 ). Because there 
is no overlap of the band with the K band and most of the kinematically 
allowed K (1430) region, there is a relatively small loss of information from 
the anti-0 cut. The technique used to correct for residual <ft contamination is 
described in Section 5. It is important to establish to what extent undiscovered 
contributions in the K + K~ channel such as /o(980)7r + could influence these 
result 

sODFor example, the / (980) 7r + contribution was found [3] to be a major 
contributing channel in the related Df — > K + K~tt + decay. Uncertainty in the 
amplitudes describing the K~K + channel will form the major systematic error 
in this analysis. Throughout this paper, unless explicitly stated otherwise, the 
charge conjugate state is implied when a decay mode of a specific charge is 
stated. 



2 Evidence for a small ao(980) contribution based on an analysis of FOCUS data 
was reported in a conference proceeding [4]. 
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2 Experimental and analysis details 



The data for this paper were collected in the Wideband photoproduction ex- 
periment FOCUS during the Fermilab 1996-1997 fixed-target run. In FOCUS, 
a forward multi-particle spectrometer is used to measure the interactions of 
high energy photons in a segmented BeO target. The FOCUS detector is a 
large aperture, fixed-target spectrometer with excellent track resolution and 
particle identification. Most of the FOCUS experiment and analysis techniques 
have been described previously [7,8,9,10]. To obtain the signal displayed in fig- 
ure 1, we required that the K~K + tt + tracks formed a vertex with a confidence 
level in excess of 10%, and the K~K + tt + vertex was outside of our BeO target 
and all other spectrometer material by at least 3a. Other tracks along with 
information from the K~K + ir + "secondary" vertex were used for form a "pri- 
mary" vertex and the separation between the primary and secondary vertex 
was required to exceed 7a. The kaon hypothesis was favored over the pion 
hypothesis by 3 units of log likelihood in our Cerenkov system for the K~ and 
K + candidates, and the 7r + candidate response was consistent with that for 
a pion. Finally, we required that the secondary vertex was well isolated. We 
required that no primary vertex track was consistent with the secondary ver- 
tex with a confidence level exceeding 0.5 % , and no other track in the event 
was consistent with the secondary vertex with a confidence level exceeding 
1 x 10~ 4 . In this analysis we eliminated most of the D + — > </>7r + contribution 
by imposing an anti-0 cut that required ttik+k- > 1050 MeV/c 2 . 
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Fig. 1. The K~K + ir + invariant mass spectra for all of the cuts used in this analy- 
sis. The D + — » K~ K + n + signal and sideband regions are shown with vertical lines. 
The peak to the right of the D + -> K'K + ir + peak is due to D+ -» K~K + tt + . 
The Df peak appears over a broad, ramped background due to misidentified 
D + — > K~ir + ir + decays. 

We obtained a signal of 6400 D + — ► K~K + tt + events prior to the anti-0 cut 
and 4200 events after the anti-0 cut. 
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3 Projection Weighting Technique 

The projector method used in this analysis is nearly identical to that used to 
determine the q 2 dependence of the helicity form factor in D + — > K~tt + £ + v 
decay [5]. To apply the projector method, we must assume that after we impose 
the anti-0 cut, the residual effects of the K~K + amplitude contributions are 
small enough that we can reliably correct for them using the bias technique 
described in Section 5. In the absence of K~ K + resonances, we can write the 
decay amplitude in terms of mx- ff + = m and the decay angle 9 (which is the 
angle between the K~ and the K + in the K~n + rest frame). 

s,p,d--- 

A = A e (m) d e 00 (cos6) 

e 

Here <iQ (cos#) are the Wigner d- matrices that describe the amplitude for 
a K~ir + system of angular momentum t to simultaneously have angular 
momentum along its (D + frame) helicity axis and the K~n + decay axis. For 
simplicity, we illustrate the technique assuming only P-wave and S-wave con- 
tributions although contributions up to and including D-wave are included in 
the analysis. In the absence of D-wave or higher amplitudes, the decay inten- 
sity consists of three terms which depend on the two complex functions: S{m) 
and P{m) of m = m K - n +. 

I(m, cos 6) = \A\ 2 = | S(m) + P{m) cos6 | 2 = 
\S(m)\ 2 + 2Re{S*(m) P(m)} cos 9 + \P(m)\ 2 cos 2 6 

For notational simplicity we will write the direct terms as SS(m) = ^(m)) 2 
and PP(m) = \P{m)\ 2 and the interference term as SP(m) = Re {S*(m) P(m)}: 

I(m, cos 6) = SS (m) + 2SP (m) cos 6 + PP{m) cos 2 6 (1) 

Our approach is to divide cos 6* into twenty evenly spaced angular bins. Let 

l D = (^m i n 2 ... *ri2o ^ be a vector whose 20 components give the population 

in data for each of the 20 cos^ bins. Here * specifies the i th m K - 7T + bin. Our 
goal is to represent the l D vector as a sum over the expected populations 
for each of the three partial waves. We will call these 'm vectors. For this 
simplified case, there are three such vectors computed for each m K - n + bin: 
{ l rh a } = ^rhss/msp/mpp). Each l m a is generated using a phase space and 
full detector simulation for D + — > K~K + ir + decay with one amplitude turned 
on, and all other amplitudes shut off. Hence the l fhpp vector is computed 
assuming an intensity of cos 2 9 for each mK-n+ bin. The 'm vectors incorporate 
the underlying angular distribution as well as all acceptance and cut efficiency 
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effects. In particular, the anti-0 cut creates substantial inefficiencies at low 
cos 6* and high m K - n +, whereas the other cuts have a reasonably uniform 
acceptance in cos#. 

We use a weighting technique to fit the bin populations in the data to the 

. — * . . 

form: l D = Fss (m{) l fhss + Fgp (m,) l rhsp + -Fpp (m{) % fhpp. The term 
Fsp(mi), for example, is proportional to SP(rrii) along with the overall ac- 
ceptance and phase space for an SP interference term generated in a given 

. — * 

rrii bin. Multiplying the % D data vector by each ! m„ produces a "component" 
equation: 
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The formal solution to this is: 



Fss (m^ 
F S p (m^ 
F PP (rrii) 



\ ( - - - - - - \ 1 ( i 

m S s ■ m ss m ss • m S p m S s • m PP 1 



m S p ■ m ss m S p ■ m S p m S p ■ m PP 
rhpp ■ rhss rhpp ■ rhsp rhpp ■ rhpp 



rhss ■ A 
rh S p ■ D, 
y l rh PP ■ Di J 



. — * — * . — * 

This solution can be written as Fss (mi) = l Pss • A, F S p (m^ = l Psp 

. — * — * 

and Fpp (rrii) = % Ppp • A , where the projection vectors are given by: 



A 



(B \ 

^ss 

Psp 

ypppj 



m ss ■ m ss m ss • m SP m ss • m PP 
m S p ■ m ss m S p ■ rh SP rh SP ■ rh PP 
rhpp ■ rhss mpp ■ rhsp rhpp ■ fhpp 



m ss 
l rh S p 
% fhpp 



(2) 



We can modify the projector weights to take into account cut efficiency, ac- 
ceptance and phase space corrections required to convert say Fsp(mi) into 

. — * 

SP(rrii) as discussed in Reference [5]. This procedure scales the % Psp , for 
example, into a modified weight l P' S p which includes all acceptance, efficiency 
corrections and phase space effects. 

The various projector dot products are implemented through a weighting tech- 
nique. For example, if we are trying to extract the 2 S(rriK-n+) x P( m K--K+) 
interference in the i th m K - n + bin, we need to construct the dot product: 



i pi .in 
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We can do this by making a weighted histogram of rriK-n+ where all events 
that are reconstructed in the first cos 9 bin are weighted by * P' SP 1 ; all events 

that are reconstructed in the second cos 9 bin are weighted by 



P' 



etc. 



4 Amplitude ambiguity 



Table 1 shows the cos 9 dependences for each of the 3 direct and 3 interference 
terms for the S(m K - n +), P(m K - 7T +), and D(rriK-n+) amplitudes. 

Table 1 

The angular distributions for amplitudes up to and including D-wave. Here cos 9 
is the angle between the K + and K~ in the K~ir + rest frame. These terms are 
products of the Wigner d- matrices d,Q (cos9). 



SS 


1 


PP 


cos 2 9 


SP 


cos 9 


PD 


cos 9 (3cos 2 (9-l)/2 


SD 


(3 cos 2 9- 1) /2 


DD 


(3cos 2 t9-l) 2 /4 



The problem is that these six terms are not independent. The relationship 
between them is given in eq. (3). 



m S D = 



3 l fhpp 



m S s 



(3) 



Hence one cannot make independent projectors of SS, PP, and SD, and must 
choose two out of these three or the inverse matrix of eq. (2) will become 
singular. Since it is known [3] that the K~n + spectrum in D + — > K~K + n + 
is dominated by the PP contribution (K ) and SS contribution [K (1430)], 
we made the choice of dropping the SD interference term and choosing SS, 
PP, DD, SP, and PD. As a result, the five projectors we use can only partially 
block any potential SD term. As long as there are no amplitude contributions 
beyond D-wave, one can show that a potential SD interference contribution 
will contaminate both the SS and PP spectra but not the two interference 
terms: SP or PD. 



The Bias Correction 



Figure 2 compares the five reconstructed spectra with the input spectra ac- 
cording to a simulation of the E687 model [3]. The discrepancies between the 
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Fig. 2. The first points ("crosses") with error bars show the reconstructed spectra 
using the method described in Section 3. The closely adjacent points ("diamonds") 
are the actual m K - n + spectra used in the simulation based on the model of Refer- 
ence [3]. These two Monte Carlo spectra are normalized near the peak bins of the 
prominent K present in the PP spectrum. The plots are: (a) S 2 (?n_ft:- 7r +) direct 
term, (b) 2 S{m K -^+) x P(m K - 7T +) interference term, (c) P 2 (m K - 7T +) direct term, 
(d) 2 P(m K - n +) x D(m K - n +) interference term and (e) D 2 (m K ~ n +) direct term. 

reconstructed and input spectra are due to residual <fi contamination in the tail 
beyond our 1050 MeV/c 2 anti-0 cut. The discrepancies are most prominent at 
the high end of the m K - n + distribution (where the <p tail is the largest), but 
are relatively small on the scale of our statistical error elsewhere. 

In order to correct for the <p tail that extends past ttik+k- — 1050 MeV/c 2 , we 
take the difference between the simulated input and reconstructed spectra as 
a "bias" that we subtract from the data after normalizing the bias by the ratio 
of the K*° peak area of the PP spectra in the data to that in the simulation^] 
Varying the K~K + amplitudes relative to the E687 model [3] will change 
the bias correction. Uncertainty in these amplitudes is our major source of 
systematic error. 



3 The normalization is actually done by finding a ratio that minimizes the x 2 
between the data and MC prediction at the peak and one bin on either side of 
the peak. 
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6 Systematic Errors 



The major source of systematic error in this analysis is due to possible uncer- 
tainty in the K~K + system amplitudes. In Section 5 we discussed the method 
used to correct the amplitude spectrum in K~n + for residual K~ K + contribu- 
tions past our anti-0 cut of itlx+k- > 1050 MeV/c 2 . This method, however, 
depends on our model for the K~K + channel. The nominal results assume 
the 4> amplitude measured in the E687 [3] analysis, but we have considered 
variations in the result due to differences in the <f) parameters as well as po- 
tential contributions from the /o(980). Although the Df — > K~K + n + has 
an 11% contribution from / (980)7r + , there was no evidence for an /o(980)"7r + 
contribution in D + — > K~K + tt + in the E687 analysis. 

We have considered possible / (980) contributions or <fi parameter variants 
consistent with our data in four areas: (a) the fraction of D + — > K~ K + n + 
that appears as <pn + : (b) the shape of the cos#£- +7r+ distribution in the vicinity 
of the 0: (c) agreement with the shape of the complete m K + K - spectrum and 
(d) agreement with the observed populations in two dimensional bins of the 
D + — > K~ K + -k + Dalitz plot. Here cos 0* K+7T+ is the cosine of the angle between 
the K + and the n + tracks in the K + K~ rest frame. In each case, the default 
model with a wide K (1430) and no /o(980) amplitude was among the best 
models matching our data according to these tests. Essentially the same set 
of potential fit variants were selected by the four tests. 

We found that a potential /o(980) when inserted with a phase of n/2 relative 
to the K (890) satisfied our consistency criteria even with a fit fraction as 
large as 3.3%. Such large / (980) amplitudes would be more inconsistent with 
the data if they came in with a relative phase of zero. We construct the 
systematic error as the r.m.s. spread of all acceptable fit variants with / (980) 
amplitudes up to 40% of the K amplitude that both matched the observed 
binned Dalitz plot populations to within 10<r of the best model and matched 
the observed cos 0* K+1T+ distribution in the vicinity of the with a confidence 
level exceeding 1 x 10~ 4 . Our default model matched the cos 0* K+n+ distribution 
with a confidence level of 18%. Table 2 in Section 7 compares the sizes of 
systematic and statistical errors for our five amplitude products. Generally 
the systematic error was found to be smaller than the statistical error. 

A variety of other checks of the results were made during the course of this 
work. We generated a Monte Carlo using the mK-w+ amplitudes summarized 
in Table 2 and compared the simulated to the observed cos 8 distributions as 
a function of m K - 7T + as well as the simulated and observed m K + K -, m K - 7T +, 
and m K + 7T + global mass projections. Agreement was good. We also compared 
the fit results obtained by analyzing the data in different ways. For example, 
we analyzed the data by constructing only three rather than five projectors 
and found consistent results with the S 2 {rriK-n+) , P 2 { m K--w+)-, and the 
2 S{rriK--K+) x P{m K - 7T +) interference terms. We raised the anti-0 cut from 
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i^k+k- > 1050 MeV/c 2 to m K + K - > 1100 MeV/c 2 and found consistent 
results even though the errors in the high mass bins went up dramatically 
unless the D-wave projectors were excluded. 



7 Results 

Figure 3 and Table 2 show the results of this analysis. The relative correlation 

between the 5 amplitudes is typically less than ±40% except at the highest 

m K - n + bins where they are as large as ±65%. The comparison plot is based 

on the E687 analysis [3] but with a much wider K (1430) represented as a 

Breit Wigner resonance with a pole at m = 1412 MeV/c 2 and a width of 

T = 500 MeV/c whereas the standard PDG [1] K (1430) parameters are 

mo = 1414 MeV/c 2 and T = 290 MeV/c 2 . In order to roughly reproduce 

the results in figure 3 in plot (a), we also increased the magnitude of the 

K (1430) amplitude by 40% relative to that obtained from the E687 analysis. 

Hence the curves in fig. 3 are drawn using a model with a 41.9% K (1430) fit 

fraction compared to the 37 ± 3.5% fit fraction quoted in Reference [3]. The 
^o 

use of a wider K (1430) with a larger amplitude increased the level of the SP 
interference in the model by 45% to a level in approximate agreement with 
our data as shown in plot(b). 



8 Summary and Discussion 

We presented a non-parametric amplitude analysis of the K~n + system in 

D + — > K~K + n + decay using the technique described in Reference [5]. There is 

no need to assume specific Breit- Wigner resonances, forms for mass dependent 

widths, hadronic form factors, or Zemach momentum factors. As described in 

Section 3, a set of five weights were generated using Monte Carlo simulations 

that were designed to project the various amplitude contributions. The five 

amplitude contributions appear in just five weighted histograms in the rriK--K+ 

mass. Because this is an essentially one dimensional technique, we chose the 

D + — > K~K + ir + final state as a test case. According to an older, traditional 

Dalitz analysis done in E687 [3] the D + — > K~K +r n + final state is particu- 

larly simple consisting of just 07r + , K K + , and K (1430)i^ + . Hence only one 

narrow resonance, the 0, should contribute to the K~K + channel that can be 

significantly reduced through an anti-0 cut such as m K + K - > 1050 MeV/c 2 . 

This leaves an amplitude that depends primarily on the K~n + mass. Another 

attractive feature of the D + — > K~K + n + final state is that there should be 

a strong S-wave component in the K~ K + channel that E687 [3] modeled as 
^o 

the K (1430). Not that much is known about including broad S-wave reso- 
nances in charm Dalitz analyzes, and we indeed find considerable discrepan- 
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Fig. 3. This figure shows the amplitudes measured in this analysis including system- 
atic errors in 73 MeV/c 2 bins. The default model described in the text are the curves. 
The data has been normalized such that the P 2 {mx- 7r +) value at the K* mass is 1 
in plot (d). The plots are: (a) S 2 {m K ^ 7! +) direct term, (b) 2 S{m K - n +) x P(m K - 7T +) 
interference term, (c) P 2 {m K -^+) direct term, (d) 2 P(m K - n +) x D(m K - n +) in- 
terference term and (e) D 2 {m^- 1T +) direct term. 



cies between our non-parametric description of the S-wave, K ir + amplitude 
in D + — > K~K + n + and the standard, PDG Kq°(1430) parameterization used 
by E687 [3]. 

The other possibly significant difference in these results compared to the E687 
analysis involves the "glitch" in the first three bins of the S 2 (mK-n+) spec- 
trum. We have observed this effect in all variants of the fit we made in the 
course of understanding systematic errors. The same bins are also present as 
a glitch in the 2 P(m^- 7r +) x D(rriK-TT+) interference term which is otherwise 
consistent with zero. If these are deemed to be significant, one explanation 
would be the presence of a small D-wave component at low K~tc + masses 
which could account for both glitches through the mechanism discussed in 
Section 4 while being small enough to escape notice in the direct D 2 (m K - 7r +) 
terms where it is second order in the amplitude. Section 4 argues there is an 
ambiguity in the cos 9 distributions if one includes waves up to and including 
D-wave which cannot be eliminated by this or any other method. Hence an 
2 S{m K - n +) x D{rriK--K+) interference term will also affect the S 2 (rriK-TT+) 
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Table 2 

Results for the five amplitude contributions. The first error bar is statistical and 
the second is systematic. All data are arbitrarily scaled by a common factor such 



that the peak P 2 (m K - n +) value near the K*° peak is set to 1. 



(GeV/c 2 ) 


SS 


SP 


PP 


0.670 





038 


± 





015 


± 





002 


0.007 ± 0.020 ± 0.002 


-0.001 ± 0.037 ± 0.004 


0.743 





018 


± 





008 


± 





003 


0.018 ± 0.011 ± 0.002 


0.035 


± 


0.020 


± 


0.002 


0.816 





018 


± 





009 


± 





003 


-0.112 ± 0.014 ± 0.005 


0.198 


± 


0.025 


± 


0.007 


0.889 





048 


± 





016 


± 





006 


-0.389 ± 0.025 ± 0.021 


1.000 


± 


0.044 


± 


0.000 


0.962 





058 


± 





009 


± 





004 


-0.097 ± 0.014 ± 0.011 


0.181 


± 


0.024 


± 


0.011 


1.035 





067 


± 





008 


± 





005 


-0.043 ± 0.012 ± 0.009 


0.033 


± 


0.021 


± 


0.012 


1.108 





083 


± 





009 


± 





007 


-0.021 ± 0.012 ± 0.013 


0.011 


± 


0.022 


± 


0.015 


1.181 





112 


± 





013 


± 





011 


-0.008 ± 0.014 ± 0.012 


0.001 


± 


0.033 


± 


0.011 


1.254 





125 


± 





021 


± 





011 


-0.053 ± 0.023 ± 0.011 


0.201 


± 


0.067 


± 


0.023 


1.327 





111 


± 





029 


± 





020 


0.037 ± 0.032 ± 0.021 


0.237 


± 


0.082 


± 


0.041 







m K-TT + 

(GeV/c 2 ) 



0.670 
0.743 
0.816 
0.889 
0.962 
1.035 
1.108 
1.181 
1.254 
1.327 



PD 



0.059 ± 0.053 
-0.078 ± 0.027 
-0.107 ± 0.034 
-0.058 ± 0.061 
0.001 ± 0.031 
-0.039 ± 0.026 
-0.004 ± 0.033 
-0.031 ± 0.053 
-0.128 ± 0.099 
-0.045 ± 0.126 



± 0.003 
± 0.012 
± 0.004 
± 0.016 
± 0.017 
± 0.012 
± 0.009 
± 0.016 
± 0.022 
± 0.063 



DD 



0.034 ± 
-0.043 ± 
-0.043 ± 
-0.011 ± 
-0.011 ± 
0.017 ± 
0.063 ± 
0.011 ± 
0.015 ± 
-0.090 ± 



0.066 ± 0.007 
0.037 ± 0.004 
0.046 ± 0.021 
0.077 ± 0.027 
0.041 ± 0.025 
0.034 ± 0.009 
0.040 ± 0.007 
0.059 ± 0.006 
0.095 ± 0.018 
0.127 ± 0.028 



spectrum. 



Although the D + — > K K + tt + decay is an ideal case for an analysis of 
this kind, it might be possible to extend the analysis to the related — > 
K-K+K+ decajQas well as D° -> K + K~K°. The emphasis in D° -> K+K~K° 
decays would be on studying the itik+k- spectrum after applying cuts to min- 
imize the effects of K ± K° contributions such as the af (980). One could also 
use this technique to study the dipion amplitudes present in hadronic four- 
body decays such as D — > K + K^n + n^ — > (j)it + -K~ decay. In particular, one 
could compare the dipion spectra produced against longitudinally and trans- 
versely polarized 0's using techniques similar to those described in Reference 
[5]. 



4 Care would need to be taken in eliminating or compensating for the D + — > 
K~tt + tt + reflection shown in figure 1, as well as the /o(980) and /o(1710) contribu- 
tions [3]. 
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